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Motivated by recent experiments, we study four-wave- mixing in an atomic double-A system driven 
by a far-detuned pump. Using the Heisenberg-Langevin formalism, and based on the microscopic 
properties of the medium, we calculate the classical and quantum properties of seed and conju- 
gate beams beyond the linear amplifier approximation. A continuous variable approach gives us 
access to relative-intensity noise spectra that can be directly compared to experiments. Restricting 
ourselves to the cold-atom regime, we predict the generation of quantum-correlated beams with a 
relative-intensity noise spectrum well below the standard quantum limit (down to -6 dB). Moreover 
entanglement between seed and conjugate beams measured by an inseparability down to 0.25 is 
expected. This work opens the way to the generation of entangled beams by four-wave mixing in a 
cold atomic sample. 
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I. INTRODUCTION 



Four-wave mixing (FWM) has been identified early on 
as a very efficient process to generate intense, non classi- 
cal beams [l|]. Indeed, the first experimental demonstra- 
tion of squeezed light was made usingFWM in a sodium 
atomic beam almost 25 years ago [2|- More recently, 
FWM in three-level A systems was predicted to gener- 
ate squeezing by using a counter-propagating geometry 
[j-[g. In addition Electromagnetically Induced Trans- 
parency (FIT) in double-A atomic systems adds flexibil- 
ity to the control of the FWM process, leading to large 
parametric gain and oscillations [3, Q • With such a sys- 
tem, van der Wal et al. generated intensity correlations 
between light pulses with a corresponding 0.2 dB of noise 
reduction below the standard quantum limit (SQL) Q. 
FWM in a double- A-system in a hot rubidium vapour and 
with a co- prop agating laser geometry was implemented 
in refs. J10l - ll5l | and by some of us in ref. [lal to generate 
bright correlated beams with a high degree of intensity 
correlations. 

In order to produce paired photons / twin beams tuned 
near an atomic resonance, efficient FWM in atomic sam- 
ples with a low parametric gain requires large optical 
depths. On the other side cold atom samples provide a 
precise control of atomic parameters as inhomogeneities 
and decoherence processes. Thus, the combination of 
double-A FIT and cold atoms has allowed the observa- 
tion of FWM at ultra-low optical powers. An optically 
thick cesium MOT allowed Kuzmich et al. [13] to gener- 
ate nonclassical photon pairs with a programmable de- 
lay. Generation and control of narrow-band paired pho- 
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tons were demonstrated in a rubidium MOT Il8|j2l| by 
reaching a high optical depth (OD), around 60 [22|. Cor- 
related photons were also produced by using cold atoms 
inside a cavity |23j . 

Following [2J], recent publications (25l-l27l| have the- 
oretically explored the time correlations between pho- 
tons observed in cold atom experiments. Propagation 
of the generated photons through the optically thick 
medium and Langevin forces represent the key elements 
of these theoretical analyses. Following the seminal 
papers [3, lH, l2J| and the analysis of [2al, we investi- 
gate the generation of quantum correlated beams in co- 
propagating beam geometry. This work presents a de- 
tailed theoretical treatment of FWM in an optically thick 
medium composed by cold atoms described by a four level 
model, in a double-A configuration. 

The main target is to explore the production of cor- 
related beams in the atom-laser configuration experi- 
mentally investigated in refs. Il0l - ll6l |. but using a cold 
atom sample as in refs. [19l - l22| . While the theoreti- 
cal descriptions |25l - [27| . concentrated on the few pho- 
ton regime, the experiments discussed in refs. [I0l4l6l | re- 
quire to tackle the large photon number regime where 
a description based on continuous variables is required 
[28l |29| . Within this approach, we derive the quantum 
noise frequency spectra and quantify the quantum entan- 
glement between the generated beams. 

Several additional original points of our analysis should 
be listed. First, we investigate the pump-seed co- 
propagating geometry, a setup which enabled to produce 
large intensity quantum correlations in the experiments 
[10l - ll6| . Furthermore, in contrast to the phenomenolog- 
ical description of ref. J12l | , gain and propagation losses 
are intrinsically included in our microscopic approach. 
As the most important result, our analytical and numer- 
ical analysis predicts that a very large degree of intensity 



quantum correlations can be achieved also in a cold atom 
FWM experiments, comparable to that achieved in refs 
[lOl [ill . Ilq . We also investigate the ability of such a 
medium to simultaneously generate large phase anticor- 
relations, leading to entanglement, a key ressource for 
quantum information [29|. 

The paper is organized as follows. Section |ll] intro- 
duces the microscopic model as described by Heisenberg- 
Langevin equations. Section IIIII discusses the analytical 
solution for the propagating quantum fields. In Section 
II VI we derive the mean values of the generated beams. We 
introduce and calculate in Section |V] the intensity quan- 
tum correlations and phase quantum anti-correlations 
spectra between seed and conjugate beams including the 
Langevin forces contribution. Finally, in Section IVH we 
discuss the effect of decoherence processes and compare 
the results to hot vapor experiments. 



II. MICROSCOPIC MODEL - 
HEISENBERG-LANGEVIN EQUATIONS 




FIG. 1: (color online) Level diagram with relevant detun- 
ings. The two large red arrows indicate the pump coupling 
between levels |1) and |3) (resp. |2) and |4)), with equal Rabi 
frequency SI. S represents the two-photon detuning on both 
A subsystems. The strong pump and weak seed configuration 
leads to optical pumping of the atoms into level [2). 



We consider a collection of atoms having the double 
A type 4-level scheme of Fig. [TJ The four atomic lev- 
els \u), are defined by their energies £■„, and we in- 
troduce the angular frequencies a;„„ = {Eu — Ey)/h, 
u,v £ {1, 2, 3, 4}, u> V. The populations of levels |3) and 
1 4) decay with a rate F. The atomic coherence between 
levels |1) and |2) decays with rate 7 and we neglect |1) 
and 1 2) population decays since we consider cold atomic 
samples. 

The atoms interact with three optical fields, pump 
(angular frequency ujp) and seed (angular frequency Ua) 
driving the first A subsystem (|1) — )- |3) — > |2)), conju- 
gate (angular frequency Wf,) and pump driving the sec- 
ond A subsystem (|I) — )• |4) — >• |2)). The detuning 
of the pump field from the |I) — > |3) transition is de- 
noted A. The pump and the seed are near-resonant with 
the |1) —7- |3) — > |2) Raman transition with a two-photon 
detuning 6. The double-A is closed by the conjugate and 
pump beams driving the |1) — > |4) — > |2) Raman transi- 
tion, assuming energy matching. As shown in FigU] the 
pump laser is detuned hy A + ujq + 5 from the |2) — > |4) 
transition, where wq = W21 — UI43. 

The interactions of the pump beam with the atoms are 
described semi-classically. For the sake of simplicity we 
limit our study to a ID model with propagation along the 
z-axis. The seed and conjugate beams are described by 
two slowly varying quantum-mechanical operators; the 
annihilation (resp. creation) operators being denoted a 
(resp. at) and b (resp. b'^). 

Following the approach of ref. [2J| , the properties of the 
medium are described by collective, slowly-varying oper- 
ators auv iz,t) averaged over small layers denoted by their 
position z. For each layer containing a large number of 
atoms Nz, we define the following population/coherence 



operators a^i 
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where 0)31 = 0)42 = LOp, 0)32 = Wa, W41 = Wb, 
0)21 = <^43 — ^p — i^a ■ The fc„„ are the projec- 
tions of the wavevectors kuv on the z axis with ^31 and 
^42 equal to the pump wavevcctor, ^32 equal to the seed 
wavevector, and fc4i equal to the conjugate wavevcctor, 
^21 = ^23 ~ ^13 £iiid fc43 = fc42 — fc32, with the convention 
h — —h 

We write the laser-atom interaction Hamiltonian V in 
the rotating wave approximation, omitting for simplicity 
the z and t dependencies of the atomic operators: 



hN f 
V = — — / dz[(a;o-l- A + 5) (744 + A(733 + (5(T22 (2) 
^ Jo 

+gaa{z,t)a32 +gbb{z,t)a4i + -^(0-31 +0-42)] + H.c, 

where N is the number of atoms in the quantization vol- 
ume V of length L. ga.t are the atom - field coupling 
constants: ga ~ ^^jr^ and gi, = BiiSt where puv is the 
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is the clcc- 



transition dipole moment and Sa,b ~ 
trie field of a single photon. In the following, we consider 
ga = gb = g- Finally, il ~ -j— is the pump laser Rabi 
frequency with E the pump laser electric field amplitude. 
We assume the transition dipole moment p to be equal for 
the two pump-driven transitions |1) — >■ |3) and |2) — )■ |4). 
As described in ref. [2^ , the evolution of atomic opera- 
tors is determined by the following Heisenberg-Langcvin 



equations: 



d \ i ~ 

TT + luv j (Juv = -r [V, auv\ + Tuv + Kv , (3) 



where juv are the atomic dephasing rates and f^v are 
the source terms produced by spontaneous emission as 
defined in [2g|. Following the ref. [25|, the Langcvin 
operators F^v are characterized by 



Fuv{z,t))^0, 



(4) 



and 



(^Fl{z,t)F^,^,{z\t')'j = 2D,,,^^,,,S{t-t')S{z- z'). (5) 

Equation ^ defines the 256 diffusion coefficients 
Duv,u'v' ', we show in Appendix IB] that only a few of them 
will be required in next section. 

Following the approach of refs. [a, [3, [2a], the descrip- 
tion of the atom-laser system is completed by a set of 
nonlinear, coupled differential equations describing the 
propagation and temporal evolution of the quantum field 
operators 



d d\ 

^+c^ja(2,t) = igAfa23iz,t), (6) 

^+c^jSt(z,i) = -igMdu{z,t), (7) 



where M = N/V is the atomic density. 



III. ANALYTICAL SOLUTION 

The system evolution is described by a set of non- 
linear, coupled differential equations for quantum fields 
and atomic operators which require further approxima- 
tions to be solvable analytically. We assume a pump 
beam highly saturating the medium and very intense 
compared to the seed and conjugate beams. Within these 
assumptions the populations cth, (T22, ctss, (744 and the 
coherences (T31, (T13, (T42, (724 are dominantly driven by 
the pump beam. In this case the complete non-linear 
differential equations system can be separated into two 
subsystems. In the first subsystem (Eq. A.l), the effect 
of seed and conjugate can be neglected and these equa- 
tions are solved in the steady state. The correspond- 
ing solutions are given in Appendix [X] These solutions 
are then injected in Eq.(Bl) which determine the coher- 
ences 0-14, (741, (723, (732,1712, (721, <734, (743 as a function 
of a and b^ as described in detail in Appendix [B] 

Eqs. ^ and ([7]) can then be written in a closed 
form and solved analytically in the Fourier space. The 
Fourier transforms of a(i, z), 6^(t, z) will be denoted 
d{uj,z), 6^(a;,z), where for simplicity we have used the 
same notation for operators in the time and frequency 



domains. Notice our notation for the Fourier transform 
of an operator and its conjugate 

/>oo />oo 

a{uj) = / a{t) e*"* dt ; a^uj) = / at(i) ^^-^t ^^^ ^g^ 



leading to [a(a;)]^ = d^(— w). The Fourier transforms of 
the quantum operators satisfy then the following equa- 
tions: 



J_d_ 
fc32 dz 

1 d 



a{uj) = TJaiL0)a{Lu)+Ka{L0p{0j)+f^'-iL0), (9) 



——b\u) = KbHaH + 776HS^H + ^6t*-H,(io) 

where rjaiio) (resp. rji,{uj)) is the complex refractive in- 
dex of the medium " dressed" by the pump laser for the 
mode d (resp. w). Ka{uj) (resp. K),{uj)) is the parametric 
conversion coefficient from mode d to w (resp. b^ to a), 
J^^*' and T^f' are the Langevin terms originating from 
the atomic Langevin forces and can be calculated from 
Eq. (IbToI) . 

Following a standard approach of quantum optics [3C| , 
we use the input/output formalism to compute the quan- 
tum operators of seed and conjugate fields. We introduce 
the input operators, dm and 6|„ for z — and the output 
operators dout and b]^^^ for z = L. The formal solution 
of the propagation equations (Pl- pU)) is given by 



dout{uj) 



'A{lu) B{uj)' 
C{uj) D{uj) 



bl 






Fa,out{<Jj) 
Pb,out{^) 



(11) 

where the coefficients A{u!),B{uj),C{u!),D{uj) and the 
Langevin terms Fi)_out{^),Fa,out{'^) are defined in Ap- 
pendix [C] Eg. (ITT]) allows us to compute both the output 
fields mean values {(dout) , (bout)) and their fluctuations 
iSdout,Sbl^t). 

As in typical FWM experiments, in the following we 
will suppose that the input on the conjugate mode b is the 
vacuum ((&„) — 0) and that the input on the seed mode 
is an arbitrary field of mean value (din)- The average 
values of the output seed and conjugate fields are given 

by 



(ao„t(w = 0)) = A{0){din) 
bLti^ = 0)) = C(0)(a,„) 



(12) 
(13) 



We define then Ga = |^(0)P as the seed gain and Gb = 
|C(0)p as the ratio between the conjugate output and 
seed input intensities. 



IV. CLASSICAL PROPERTIES OF THE 
MEDIUM 

The mean value expressions for the seed and conju- 
gate output operators allow us to calculate the output 
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FIG. 2: Spectra of the Ga seed beam gain in a) and of the 
conjugate beam gain Gt in b) as a function of the two photon 
detuning S. The parameters used for the plot are : ^/2n — 
10 kHz; Q/2-K = 0.3 GHz; A/2-k = 1 GHz. The seed beam dip 
for (5 ~ A in a) is originated by absorption on the |2) — >■ |3) 
transition. Expanded views of Ga and Gt gains around 5 = 
are shown in the insets. The peculiar line-shapes are the result 
of a nontrivial interplay between the two processes discussed 
in the text. 



intensities which can be experimentally measured. The 
parameters used in the following for laser intensities and 
detunings are derived from hot vapors experiments such 
as pi>jl6l |. so that in our numerical calculations, we will 
use 0.7 < A/27r < 3 GHz and 0.3 < n/2TT < 2 GHz. 
We consider the Di transition of ^^Rb: lasers wave- 
length 795 nm; ljq/27t ~ 3 GHz; natural lincwidth 
r/27r = 5.7 MHz; resonant absorption cross section 
10-9 ci 



CTO 



31 



For the optical depth defined by 
OD = McToL, we use the value of 150 reachable by im- 
proving a cold-atom configuration such that of ref. [22] . 
In cold atom experiments, magnetic field inhomogeneity 
and atomic collisions may determine the 7 decoherence 
rate. In a state-of-the-art experiment, 7 depends mainly 
on the relative phase stability of the lasers and is as low 
as 7/27r ~ 10 Hz [S^. However, as specifically discussed 
in Sec. lVI Al our results shows a weak dependence on this 
parameter and 7/27r = 10 kHz will be usually assumed. 
Fig. [5] reports numerical simulation results for the seed 
and conjugate gains as function of the two-photon de- 
tuning S. As described in ref. [6| different elementary 
processes contribute to the system dynamics. Effective 
coherent photon redistribution from pump to seed and 
conjugate occurs when the four-photon resonance con- 
dition is fulfilled (S ~ 0). Let us note that the same 



resonance condition holds for a Raman absorption pro- 
cess involving the seed beam. While it is quite obvious 
that mode b will always be amplified as soon as the non- 
linearity exceeds the absorption loss, the behaviour of 
mode a is more complicated and relies on the interfer- 
ence between redistribution and Raman processes. The 
spectrum profile of the seed gain (Fig. [5^) is character- 
ized by absorption at ^ ~ A, with a width imposed by F 
and the propagation through the optically thick medium. 
For J ~ 0, the seed gain spectrum displays a sharp pro- 
file containing both gain and absorption contributions 
(Fig. [3i inset). The seed absorption is due to the Ra- 
man process on the transition |2) — >■ |1) involving one 
seed photon and one pump photon, while the gain can 
be attributed to the FWM process. For a pump laser in 
resonance with one "arm" of the A transition the KIT 
profile assumes the characteristic lineshapc of a narrow 
dip inside a Lorentzian profile. In an KIT configuration 
with the pump laser detuned from the resonance, the 
absorption profile of the seed field becomes asymmetric 
about the two-photon resonance assuming a characteris- 
tic Fano-like profile as in the experiments of refs. [33l436l | 
and the theoretical analysis of ref. [33] . The i5 ~ pro- 
file of Ga is similar to the EIT Fano profiles reported in 
the above references, except that Ga > 1 in a narrow 
frequency range. The precise positions of the seed ab- 
sorption and Fano-like feature are given by the AC Stark 
shifts induced by the pump laser. 

The Gb spectrum, shown in Fig. [2]d reports a reduction 
of the gain in the region (5 ~ with a width matching that 
of the dip in the Ga spectrum. In a naive description, 
conjugate gain does not increase along the propagation 
when the seed is completely depleted. 



V. QUANTUM PROPERTIES OF THE SYSTEM 

A. Noise spectra 

Now we focus our attention to the quantum proper- 
ties of the system. We follow the standard approach of 
the continuous variable formalism [30] to calculate the 
noise spectra. For the seed mode a we introduce the 
following amplitude and phase quadrature fluctuations 
components, Sxa and Spa respectively: 

6xa = Sae-^'^'^ + 5a^ e^'^'^ (14) 

5pa = -i{5ae-'^'^-5a)e"f'^) (15) 

where (pa is the phase of the mean field (d) . We define ac- 
cordingly the amplitude quadrature fiuctuations for the 
conjugate mode 6, 5xb and Spb- Following the Wiener- 
Khintchine theorem, the spectrum of a given quantity u 
can be evaluated by taking the Fourier transform of its 
autocorrelation function ^Ol : 



(16) 



Su{uj)2ti5{u - J) = {u{u)u^ [uj')) 

We analyze the FWM configuration of refs. [lOl - [la ] 
with a coherent state for the seed input of mean value 



amplitude a. In agreement with the symmetric ordering 
of the operators [301 ^^ be used hereafter, the covariance 



J 



matrix for the input seed field is given by 



(5xa,,„(a;)(5x^_„(w')) {|(5.Ta,,„(a;)(5pt^,„(a;') 

I 



2'n:5{ijj + uj') 



1 
1 



(17) 



and a similar definition applies for covariance matrix of 
the conjugate mode but with a zero mean input power. 
Notice that, according to [33, 123, the choice of order- 
ing of the operators should not influence the calculated 
spectra. 



J 



Using Eq. dill) for the propagation of the quantum 
operators for the seed and conjugate fields and the input 
covariance matrix defined in Eq. (J17p , the quantum noise 
spectra for the quadrature amplitude of the fields are 
given by 






\{ao^t)\' 



[\A{^)\'{1 + i?aat (^)) + \A{^^)\'{1 + D,u{-^)) 



+ \B{uj)Wl + T^bt^M) + \B{-Lo)Wl + Z?6bt(-c^))] , 

%^ [|CHp(l + i?„,t M) + |C(-^)p(l + D,,,{-u^)) 

+ \D{lu)Wi + D^,^{uj)) + \D{-uj)Wl + D^^.i-u))] , 
I 



(18) 
(19) 



where the calculation method of the four atom-driven 
Langevin diffusion terms D„,aU -D^jta ^6tb and D(,tt is 
presented in Appendix iDl All those coefficients are, by 
construction, real and positive quantities and act as ad- 
ditional noise source terms on the quantum noise spectra. 
Their specific contribution will be discussed below. 

The interest for the double-A FWM system stems in 
particular from its ability to generate quantum field cor- 
relations. Such correlations can concern a single quadra- 
ture of seed and co njug ate fields (e.g. intensity correla- 
tions as in refs. [io, |4l[) or two conjugated quadratures 
(e.g. intensity correlations and phase anti-correlations). 
In the second case, seed and conjugate fields can be en- 
tangled Let us note that while most FWM experiments 
investigated intensity squeezing, Ref. [l4{ reports entan- 
glement in a hot atomic vapour. In order to detect the 
entanglement, we will here use the I{uj) < 1 sufficient 
criterion introduced in [42], |4J] and based on the insepa- 
rability parameter defined by 



iH = l{s- + s+). 



(20) 



Here S~ is the intensity correlation spectrum and S^ the 
phase anti-correlation spectrum both normalized to the 
standard quantum limit (SQL). S~ is derived by applying 
Eq. dHl) to 



1 

71 



iSXa,out{^) - Sxb^outi'^)) 



(21) 



and normalizing it to the sum of gain on mode a and 
mode b: (Ga + Gb)- Correspondingly 5*+ is obtained 
with 



u — 



—7= {Spa,outii^) + Spb,outiuj)) (22) 



with the same normalization. Inseparability will be used 
below to give a higher bound for entanglement as dis- 
cussed in ref. [ij . Eq. (1^0)) leads to the following expres- 
sions for the normalized intensity correlation spectrum 
and the normalized phase anti-correlation spectrum: 



S^M 



5+H 



1 



2{Ga+Gt) 



({|A{0)M(a;) - C(0)*CH|2)(1 + D„„t H) + {|A{0)A{-a;)* - C{0)C(-a;)* \^)(1 + D^^J-u;)) 



(|A(0)*BH - C(0)*DH|2)(1 + D,t,H) + {\A{0)B{-ur - C(0)D(-^)* \^)(1 + D,,, (-^))) (23) 

^,„ \„ , ((|A(0)CH - C(0)* AH* p)(l + D„„t H) + {|A{0)C{-a;)* - C(0)*A(-a;)* p)(l + D^,^{-uj)) 
(|A{0)D(a;) - C(0)*iJM*|2)(l + D^t^H) + {|A{0)D(-^)* - C(0)*B{-a;)*|2)(l + D,,^{-oo))) (24) 



In contrast to the case of Eq. p^ , which is a sum of pos- 
itive terms, a non- vanishing set of coefficients A, B, C, D 
could exist in the Eq. ((23)) . such that the spectrum 
of the intensity difference S~{uj) —>■ (same holds for 
Eq. 24). For example, in the case of a system be- 
having as an infinite bandwidth ideal linear amplifier 
having gain equal to G (|A(w)p = |D(a;)p = G and 
|B(w)p = |C(a;)|2 = G - I, and vanishing Duv), the 
relative intensity noise spectra is given by 

I 



^."M = 



2G-I' 



(25) 



as reported in [l2| • Let us note that our model describes 
the frequency dependence of the noise spectra as well as 
the Langevin forces contribution through the coefficients 

DaaU -Data i'fctfo and D^bi . 

Fig. [3] displays the intensity correlations and phase 
anti-correlations spectra versus a;/27r (analysis fre- 
quency) for two values of O and A. For each set of values 
the two-photon detuning S is chosen in order to optimize 
the noise reduction, close to the maximum gain value. In 
both cases the output beams are entangled since they dis- 
play both non classical intensity correlations and phase 
anti-correlations. Our calculation shows that, within the 
explored range of parameters, at fixed fl/A the optimum 
entanglement increases with fl. Let us note that both 
seed and conjugate output beams display an intensity 
noise above the standard quantum limit as reported in 
hot atomic vapor experiments [lO|. 



B. Contribution of Langevin forces 



Eqs. ([23)) and p4|) show that the contribution of the 
Langevin forces is purely detrimental for entanglement. 
In a previous paper studying similar systems |26| . the 
contribution of Langevin forces was negligible. In our 
system we found that it was not the case and we have 
specifically investigated their role. Fig. 2] shows the in- 
separability spectra versus the analysis frequency in pres- 
ence and absence of the Langevin terms. Their effect is 
small but clearly non- negligible and increases with w. 
Let us point out that, when neglecting the Langevin cor- 
rection terms, our calculation leads to non-physical spec- 
tra for the mode a, i.e., the noise can fall simultaneously 
below the SQL for both phase and intensity quadratures. 
On the contrary, by including the Langevin terms, we 
always obtain, as mentioned above, excess noise on both 
quadratures. 



VI. DISCUSSION 

Below, we discuss briefly the role of the ground state 
decoherence rate as well as the comparison of the model 
to the results of a hot vapor experiment. 
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FIG. 3: (color online) Intensity correlations (black dashed) 
and phase anti-correlations (red continuous) spectra (in dB) 
as a function of analysis frequency for a) Q./2ti — A/27r = 
0.3 GHz and S/2tt =-48 MHz and b) n/27r = A/27r = 2 GHz 
and S/2tt =-217 Mhz. Inseparability can be obtained by tak- 
ing the half sum of the two correlation spectra. 
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FIG. 4: (color online) Inseparabilty with (red line) and with- 
out (black dashed) the Langevin noise correction. Parameters 
as in Fig. 3b). 



A. Role of the ground state decoherence rate 

Quantum correlations are notably sensitive to various 
decoherence mechanisms. In Fig. [5] we present the evolu- 
tion of the maximal quantum correlations (in terms of the 
I inseparability) as a function of the 7 decohernece rate. 



While high relaxation rates (7/27r > 10 MHz) do not al- 
low for the observation of entanglement, the smooth be- 
havior in the low 7 region shows that the system is quite 
robust against mild decoherence mechanisms. These ob- 
servations show that coherences between hyperfine levels 
play a key role for the production of entanglement in 
FWM experiments. 



by the cold-atom simulation. 
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FIG. 5: Effect of the decoherence between the two ground 
state on the Tiyi) inseparability at u)/2ti = 1 MHz. One can 
see that the process is robust to decoherence up to 1 MHz. 
Parameters: n/27r = A/27r = 2 GHz, 5/27r = -217 MHz. 



B. Comparison to hot atomic vapor experiments 
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FIG. 6: Gain spectrum Ga of a weak seed beam as a function 
of the two photon detuning 5. The black dashed line reports 
the experimental data for a heated Rb cell, as published in 
ref [U by the NIST group (courtesy of P. D. Lett). These 
results are obtained with a pump power P=400 mW, a pump 
beam waist of 650 fj,ra {Q,/2-k — 330 MHz), a pump detuning 
A/27r = 700 MHz, and L = 12.5 mm. The red solid line 
corresponds to the theoretical predictions of our model using 
those parameters and an atomic density A/" = 4 x 10^ cm^ in 
agreement with that estimated from the measured rubidium 
cell temperature. The gray shaded area corresponds to the 
spectrum reported in Fig. [2^ (please note the reversed scale). 
The absorption dip and the Fano-like profile on the right part 
of the figure, absent in the simulation of Fig. O are originated 
by the seed acting on the |1) — !> |4) transition. 



The presented numerical results were obtained assum- 
ing a stationary regime (no transient phenomena asso- 
ciated to transit-time) and in the absence of inhomoge- 
neous Doppler broadening. The extension of the model 
in order to include these effects is beyond the scope of 
the present work. However, let us present the result of 
the calculation of the Ga seed gain, using our cold atom 
model, assuming the set of atom/laser parameters ex- 
plored by McCormick et al [10|. For the atomic density, 
we assume Af ~ 4.10^^ cm~'^ falling into the range ex- 
plored by the experiments of refs. (l(]| - ll6l |. We estimate 
the decoherence rate 7/27r — 500 kHz from the mean 
transit time through the pump beam. 

Fig. ini shows our numerical result for Ga spectrum and 
compares them to the experimental observations kindly 
provided by P. D. Lett. Notice that the gain spectrum 
contains twice the absorption and Fano-like feature of 
Fig. 2, because the seed and conjugate beams exchange 
their roles in a large scan of the seed laser frequency. 
Even if the cold-atom model cannot reproduce exactly 
the data obtained in a heated cell where the Doppler 
effect plays a major role, the main features of the exper- 
imental transmission profile are well reproduced. In ad- 
dition the gain peaks occur at the predicted two-photon 
detuning; the coherent absorption dips are present and 
finally the asymmetric profile of the left peak is also ac- 
counted for. The gain peak values of the hot vapor exper- 
iment are in very good agreement with those predicted 



VII. CONCLUSION 

We explored analytically and numerically the FWM 
generation process for an atomic doublc-A system. The 
Heisenberg-Langevin formalism allowed us to describe 
the classical and quantum properties of seed and conju- 
gate beams beyond the linear noiseless amplifier approx- 
imation. A continuous variable approach gave us access 
to relative-intensity noise spectra (that can be directly 
compared to experiments). In the very wide parameter 
space of this system, we explored a range of large pump 
Rabi frequency and detuning, recently addressed by hot 
Rb vapors experiments. However, Doppler broadening 
was not taken into account in our analysis and the op- 
tical depth was assumed compatible with a cold-atom 
sample. In this previously unexplored regime, we pre- 
dict the generation of quantum-correlated beams with a 
relative-intensity noise spectrum well below the standard 
quantum limit, down to -6 dB. Moreover, the calculations 
predict entanglement between seed and conjugate beams 
characterized by a inseparability down to I = 0.25. We 
verified that in the explored regime the Langevin forces 
cannot be neglected and reduce the performances of the 
system. We examined the role played by the ground- 
state (hyperfine) coherence in the generation of quantum 
correlations: the smooth behavior in the low 7 region let 



us predict that the system is robust against mild deco- 
herence mechanisms, as magnetic field inhomogeneities. 
Our model investigating the double A scheme is quite 
general. It may be applied to describe a wide range of 
systems having the same level structure, for instance in 
solid state physics. 
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APPENDIX 



Appendix A: Steady state 



In zeroth-order perturbation expansion, in which a 
and b^ go to zero, the Heisenberg-Langevin equations for 
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The steady state solution of Eq. (jAll) is 

|(So)] = [Mo]-i|5o]. 
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Appendix B: Atomic Heisenberg-Langevin equations 

The first order solution for the four coherences, 
<723j 1741, (743, (721 IS determined by the following matricial 
equation 



^[l]^^ + [^h]]\J:,]^\S,]\A] + ^\F,] (Bl) 
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A similar set of equations fields for the hermitian conju- 
gate operators. 

The Langevin atomic forces \F] are characterized by 
their diffusion coefficients matrix [D] (Eq. ([5])). In order 
to comply with the symmetrical ordering condition intro- 
duced in Section |Vl wc write the following symmetrize 
diffusion coefficients 



(|Fi(z,t)][Ft(z,i')l), (B5) 



[D] = {[D,] + [D,]), 



(B4) 



[D,]2S{t-t')S{z^z')^{\Fl{z,mF,{z,t')\). (B6) 



Langevin diffusion coefficients for operators can be 
calculated using the generalized Einstein relation as in 
ref. [301 • The [Di] and [D9] diffusion matrices are given 
by 
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The system of Eq. (jBTj) is solved by writing each 
atomic operator as the sum of its mean value and a quan- 
tum fluctuation term 

|I]i] = |(I]i)] + |<5S]i], (B8) 

By linearizing Eq. (jBip wc derive for the mean values 

|(Si)] = [A/i]-i|5i]|(i)]. (B9) 

and for the Fourier transformed quantum fluctuations 



|(5Si] = ([A/i]+< 



\Si]\6A]+i{[Mi]+Lo[l])-^\F, 



(BIO) 



where [il/(w)] is a 2x2 matrix given by 



[M(^)] = _,^[T]([Mi] +c.[I])-i|5i] (C2) 



with [T] = 
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10 



The Langevin force terms of Eq. pT|) are given by 

nl 
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Appendix C: Quantum field Heisenberg-Langevin 
equations 

The propagation matrix for the quantum held opera- 
tors of Eq. (|11|) is obtained by replacing the o'23(t, z) and 
(T4i(t, z) solutions of Ea. (jB8[) into the differential Eqs. (|6]) 
- ([7]) and solving these equations. By deflning 
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Let us point out that the conjugate operators for the 
quantum flelds are given by 
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where we have introduced \F^] 



-P32 
F12 



By defining the D^^^^\ (w) coefficient as 



Appendix D: Langevin forces contribution to noise 
spectra 



(F,HF,t(c.')> = D,,f{u;)2n <5(^ + ^') 



(D3) 



We will present this calculation for the term 
(Fa(w)FQt (w')) term, as the extension to the remaining obtain 

ones is straightforward. Substituting the expression of 
Eq. (|C6|) this term is given by 



(FaHF,t(c.')> = ([10| / e-[^^(")l^^[MfH]|i^i(z,c.)]dz 
Jo 



The Langevin forces are i5-corrclated in z so integration 
over z' gives : 



Jo 

x*[M*(-w')]e"*'*^*^""'^'^^dz)|l 0], (D4) 



(Fa(w)/;t(w')> = [1 0|( / e"[*^(")l'^^[A/F(cj)]|Fi(2,w)] within the [D] diffusion matrix given by Eq. ^. In a 

•^ similar way it is possible to derive D^i^{Ld), Di^i^t (w) and 

x[fI{Lz,uj')\ *[Af^(~a;')]e-'[*^'(-"')l'^^dz)|l 0]. (D2) Df,f,,{uj), as detailed in [H. 
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